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Abstract 

The dimer problem arose in a thermodynamic study of diatomic molecules, 
and was abstracted into one of the most basic and natural problems in both statistical 
mechanics and combinatoric mathematics. Given a rectangular lattice of volume V in d 
dimensions, the dimer problem loosely speaking is to count the number of different ways 
dimer s (dominoes) may be laid down on the lattice to completely cover it. It is known 
that the number of such coverings is roughly e XdV for some number A^. The first terms 
in the expansion of A^ have been known for about thirty years 

A, ~ i ln(2d) - l - . 
Herein we present a mathematical argument for an asymptotic expansion 

\ d ~ - ln(2d) 1 H + •■■ . 

2 v ; 2 8d 96 d 2 

with the first few terms given explicitly. 
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In a previous paper, [1] , we worked with tiles for which there was a weighting function 
assigned to the tiles. The weighting function satisfied a normalization condition. We then 
associated a quantity we called the "pressure" to the tilings of a lattice with such tiles. 
Loosely speaking, the result of the paper was that if we varied the weighting function, 
letting it get smoother and smoother (that is, more slowly varying), the corresponding 
pressures approached a limit. The limit was, of course, the pressure associated to the 
"infinitely smooth weighting function" a constant function. Herein we continue the sta- 
tistical mechanics approach one step further. We construct for any weighting function 
the perturbation expansion for the partition function that arises naturally from treating 
the difference between the actual weighting function and the infinitely smooth weighting 
function as a perturbation. 

We construct a second "partition function" associated to the first partition function, 
which we develop in a "cluster expansion". Kernels Ji, J2,.... arise in a standard way 
from the cluster expansion development. We next specialize to the dimer problem. We 
relate our two partition functions and we herein give a formal derivation of an expansion 
we hope in the future to prove asymptotic. 

We know in fact c\ = 1/8 and C2 = 5/96. The q are computable from the J's. The 
leading term | ln(2d) — | does indeed arise as the zeroth order term in the cluster 
expansion. 

The present work was inspired by, and is a natural sequel to, our previous paper [1]. 
However it is essentially independent thereof. We write the paper so it may be read 
without reference to [1], some notation is modified, but we strive to keep consistency 
between the two papers. 

We work with a d- dimensional lattice, A, on a torus. We let N be the number of 
vertices of A, (hopefully without confusion, we also let A be the set of vertices). We 
consider tiles, t, of size, or area, n. (Our tiles need not be connected, and N/n must be 
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an integer.) We have a weighting function, /, on tiles that is normalized so that 

t 

where the sum is over all tiles. A located tile, s a , is a tile placed at a particular location 
on the lattice, and is exactly a subset of the lattice of cardinality n, 

s a = {x 1 , x n } . (2) 

( A tile may be viewed as an equivalence class of subsets under translation, and placing a 
tile at a particular location corresponds to selecting an element of the equivalence class.) 
We denote the corresponding tile as s a . Equation (1) can be written as 

E /W = 1 ( 3 ) 

XQ£s a 

The sum in (3) is over subsets of A of cardinality n that contain a fixed point of the 
lattice x . This equation is easily identified with equation (3) of [1]; as naturally for any 
located tile s a we identify 

/fe)=/(z?,...X) (4) 
where the left side / is of this paper, and the right side / is as in [1]. 

We now consider a tiling, Tj, of A by tiles of size n. Tj may be identified with a set 
of N/n located tiles 

Ti = {s\, s l 2 , s]v/ n } (5) 

Since Ti is a tiling 

U s « = A ( 6 ) 

a 

and 

<H4 = 0, a ^(3. (7) 

The partition function, Z, and the quantity we call the pressure, p, are given by the 
equation 

z = eNp = e n ( § ) 
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The right side is the sum over tilings, each tiling weighted by the product of the weight 
functions of the tiles appearing in the tiling. 



The "infinitely smooth weight function" is constant on tiles. We let f be its value 
on any tile, easily calculated to be 

_ (n - 1)\(N - n)\ 



We write 



/(*«) = /o +(/(*«)" /o) 



(10) 



= fo + v(s a ) 

where we will treat v as a perturbation. Equation (8) is now 

Z = e Np = Y, II (/o + «(*«))■ (11) 

Ti s a &Ti 

We expand Z in terms of powers of v 

Z = Z + Zi + Z 2 + ... (12) 
Zi homogeneous of degree % in the v's. Singling out Z first 

^=En/o=E/o" /n as) 

and one calculates 

\ = (?)!(n!)^- (14) 
We put together (9), (13), and (14) to define 

Z Q = e*?™ (15) 

where 

1 — n 

lim f(N) = -. (16) 



(These are equations (8) and (9) of [1].) 



We turn for a moment to the dimer problem. [2] We relate the quantity \ d of the 
dimer problem to our variables (if our / is set equal a constant on dimers and zero on 
other tiles). 

e NX d = ( 2 qN/2 z (17) 

or 

\ d = l - ln(2d) + ^lnZ (18) 

the (2d) N/2 in equation (17) arises since our normalization condition (1) requires our / to 
be 7^7 on each dimer, whereas to just count tilings (with no weighting) would correspond 
to / = 1 on each dimer. Replacing Z by Z in (18) and taking the limit N — > oo one 

gets 

\ d & \ ln(2d) - \ . (19) 

Our expressions for A^ will all be in the infinite volume limit. Thus taking our zeroth 
order approximation for Z,Z — Z , we get the approximation for A^ of (19), the result 
of [2] ! In [2] there are bounds on the error of this approximation. 

Returning to the general case we factor out Z from Z in (12) 

Z = Z Z* (20) 

Z* = 1 + Zl + Z* + (21) 

Z* = Zi/Zo (22) 

There is a natural factoring of Z{ into a contribution from the factors of v in (11) which 
we call Z* and the factors of fo in (11) which we call (3(N,i)Zo so that 

Z l = P(N,z)Z Z* (23) 

and thus 

z; = p(N,i)z; (24) 

Z* is defined below, and it is a tedious calculation to show the surprising result that 
P(N,i) -> 1 as N -> oo. We let Z* be Z* with /3(N,i) replaced by 1. 

Z* = 1 + Z{ + Z* + ... (25) 
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We turn to the detailed specification of Z* 



Sl,S2,...,Sj Q— 1 

disjoint 

We are trying to put our expression for Z* into a form that can be identified with 
equation (2.5a) from [3]. Some knowledge of the structure of cluster expansions becomes 
necessary. Article [3] is a standard reference, but a knowledgeable patient friend may be 
more helpful. 

With the notation in [3] , the formal cluster expansion for Z* is given as a translation 
for equation (2.7) of [3]. 

ln ~ z * = E 4 Js ( 27 ) 

s 

J s = ^ v(s 1 )...v(s s )ip' c (s 1 ,s 2 ,...,s s ) (28) 

Sl,S2, — ,S a 

The located tiles appearing in the sum for J s are forced to overlap so that they cannot 
be divided into two disjoint sets of located tiles. ip' c ls a numerical factor depending on 
the pattern of overlaps. An illustrative computation of a sample J s is treated in the 
appendix. 

The J s for the dimer problem easily satisfy 

Ji = (29) 
and we will show in a succeeding paper that J s is of the form 

+ <30) 

with r > s/2. 

We want the asymptotic behavior of Z* as iV — > oo 

Z* = J2P(N,i)Z* (31) 



and we know from (27) the asymptotic behavior as N — > oo of Z* 

Z* = J2 Z* (32) 

i 

We argue the asymptotic behavior of (31) arises from a largest term of the form 

atH aiN 

P(N, £ , a, N) J- ■ ■ • W ■ (aiJf) ,... M , (33) 

provided all these J's are > 0. Here JjiV = (1/i!) Ji- We are finding that portion of the 
asymptotic behavior due to just Ji, J s +i- We choose the a, to maximize (33). We will 
use (via another tedious calculation) 

P(N,jN) ~ e"^) 1 * 1 - 2 *™] 
Differentiating (33) with respect to the ctj leads to the equations 

, 7 9 

In a k = In J k + - — 

da k 

or 



1 — 2Eia,- , 

m(l — 2E zai) + S z«i 



(34) 



a* = J fe e^ Q ' s ) (35) 

where (34) and (35) define the F k . Equation (35) is solved for a k as a formal power series 
in the J by iterating starting from setting the a's = on the right side of equation. 
Again, in a succeeding paper, we argue these same formulas hold even if some of the J's 
are negative. 

We substitute (35) into (33) to get 



JV j-E a x Fi + S Jt e F » + 1 2 \ i °' 2n(l-2£ i a^+Siaij 



~ e I 2 v v J. (36) 

From (34), (35), and (31) we can expand the exponent in (36) into a formal power series 
in the Jj, and get the asymptotic series 

l lB(M) _l + ! + ! + ... (37) 



and from 



11-11 - 3 (2d - 1) 

8 d' 3 ~ 12 d? ' 4 ~ ~ 64 d 3 



(3* 



get 



Cl = -, c 2 



5 

96 



(39) 



A2 was calculated exactly in [4] and [5]. Good bounds for A3 are given in [6]. We hope 
the asymptotic expansion we have developed actually is convergent, even for d = 2 and 
d = 3! We are actively working to show this. 

Note . In an earlier form of this paper we had an incorrect value for J 4 . (The calculation 
of J4 was indeed tricky.) We also incorrectly assumed that in equation (30) r was equal 
to s — 1. 
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Appendix 



We consider the dimer problem on a one dimensional lattice with N vertices. 
These vertices may then be labeled by the integers 1,2, N. The generic located tile 
may be written s^- where 1 < i < j < N, and clearly 

Sij = {i,j} (Al) 

/o, here, is exactly the constant function on tiles equal l/(A — 1). f is given as 



f(Sij) 



< \ if j = i + 1 

\ if i = 1 and j = N . (A2) 
otherwise 



And by (10) 

v ( s ij) = f( s ij) ~ fo (A3) 
We study J 2 = 2NJ 2 . We will actually be after J 2 in the — > oo limit. 

We refer to equation (28). We find from [3] that ip' c (s\, s 2 ) = —1 if sifls 2 is non-empty, 
and zero otherwise. 

J ^ = -\^ E' (/M-/o)(/M-/o) (AA) 

where the prime on the summation sign is to represent the restriction to s^-, Sjy with 
non-zero intersection. We write J 2 as the sum of six terms 

J 2 = A + B + C + D + E + F (A5) 

A,B, and C are terms where s^- and s^y intersect in two vertices (i.e. s^- = s^y), and 
D, E, and F are terms in which s^- and s^y intersect in a single vertex. We expand the 
binomial in the sum of (A4), and A and D are terms quadratic in the /'s, B and E, 
quadratic in the /o' s ) an d C and F are crossterms involving one fo, and one /. We want 
the limits of these quantities as N — > oo, which we here collect. 
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B — ► (A7) 

C — > (A8) 

D — -i (A9) 

E — -\ (A10) 

F — >1 (All) 

From this will follow that J 2 — > |, the result in (38), for d — 1. We derive some of the 
limits in (A6) -(All). 

Study of term B 



ij_ w-i))m- _ (412) 

The last factor on the left side of the limit is /g, and the next to last factor is the total 
number of possible located tiles. We are here keeping terms in (A4) with Sij = Si>j' arising 
from the product of the / 's- 

Study of term A 



A = -- —.N-- = -- (413) 
2 N 4 8 K J 

The factor \ is f 2 on each dimer, and the factor N is the total number of dimers. We 
are here keeping terms in (A4) with = s^f arising from the product of the f's. f is 
only non-zero on dimers. 

Study of term F 
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F = ----2-N-2-{N-2).-- w — i ^+l. {AU) 

The factor — ^ comes directly from (A4); the next factor 2 is associated with the two 
cross terms, ffo and fof] the next N arises because we pick s^- = S12 (we are looking at 
f(sij)fo, and using translation invariance of /, and the fact / is only non-zero on dimers); 
2 arises from the two choices i' = 1, i' — 2; (N — 2) is the number of choices for j', once 
i' is fixed, \ ■ -^y is the non-zero value of —ffo- 

The other terms are similarly calculable. This has been a typical computation for the 
J's. There are lots of cases; and the number of cases grows for higher J s 's, as does the 
counting become trickier for each case, as s grows. 
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